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Abstract 

Transient growth due to non-normality is investigated for the Taylor-Couette problem with 
counter-rotating cylinders as a function of aspect ratio r\ and Reynolds number Re. For all 
Re < 500, transient growth is enhanced by curvature, i.e. is greater for rj < 1 than for r\ = 1, 
the plane Couette limit. For fixed Re < 130 it is found that the greatest transient growth is 
achieved for r\ between the Taylor-Couette linear stability boundary, if it exists, and one, while 
for Re > 130 the greatest transient growth is achieved for rj on the linear stability boundary. 
Transient growth is shown to be approximately 20% higher near the linear stability boundary 
at Re = 310, r/ = 0.986 than at Re = 310, rj = 1, near the threshold observed for transition 
in plane Couette flow. The energy in the optimal inputs is primarily meridional; that in the 
optimal outputs is primarily azimuthal. Pseudospectra are calculated for two contrasting cases. 
For large curvature, rj = 0.5, the pseudospectra adhere more closely to the spectrum than in a 
narrow gap case, rj = 0.99. 



Keywords: Taylor-Couette instability, transition to turbulence 



2 



I. INTRODUCTION 



Two shear-driven flows bear the name of Couette: the flow between differentially rotating 
concentric cylinders is called cylindrical Couette flow, or more commonly, Taylor-Couette flow, 
while the flow between infinite parallel plates translating at different velocities while maintaining 
a constant separation is called plane Couette flow. Exact solutions, both called Couette solutions, 
to the Navier-Stokes equations for each of these configurations are easily calculated and serve as 
illustrations in most textbooks. Each of these two flows has become a paradigm of hydrodynamic 
stability theory. 

Taylor-Couette flow can be considered to be a paradigm of understanding. In 1923, Taylor 
|Q] carried out calculations of the linear instability of Couette flow to the onset of axisymmetric 
vortices and compared these with experiment, obtaining agreement which remains remarkable 
even by today's standards. In later research, increasingly ornate and beautiful experimental 
patterns were discovered, e.g. ||], |l| and correspondingly elaborate numerical, asymptotic, 
and theoretical calculations, e.g. ||, pi M, reproduced and explained these patterns, again with 
remarkable accuracy; see [§, ||, [l0|, [ll]] . 

Plane Couette flow, on the other hand, can be considered to be a paradigm of mystery. For 
plane parallel flows, Squire's theorem [12| firmly establishes that the linear instability with the 



lowest critical Reynolds number is spanwise invariant. Armed with this theorem, researchers 



have long known, and more recently proved [13], that plane Couette flow is linearly stable 



at all Reynolds numbers. Yet, in laboratory experiments |l4j and in numerical simulations 



[15 1, plane Couette flow undergoes sudden transition to three-dimensional turbulence. Since at 
least the 1960's researchers have explored various mechanisms for transition in plane Couette 
flow which bypass linear instability and Squire's theorem. An element shared by all of these 
approaches, independent of the theoretical mechanism proposed for instability, is the presence 
of streamwise vortices fl6|, 17, 18, [H| ^(|, i.e. perturbations which, unlike those deemed critical 



by Squire's theorem, are not spanwise invariant, but are instead invariant or almost invariant in 
the streamwise direction. These are analogous to the Taylor vortices which, in Taylor-Couette 
flow, are the eigenvectors responsible for the linear instability and are realized in experiments 
and in nonlinear numerical simulations. 

One major line of research has focused on the effect of the non- normality of the operator 
governing the linear stability of plane Couette flow jlT], 15, Such operators may lead to 
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transient growth dynamics, even in the absence of linear instability, allowing nonlinear effects 
to take over before the final exponential decay exhibited by the linear evolution. Plane Couette 
flow can exhibit transient growth of several orders of magnitude; the initial conditions which 
maximize the transient growth, called optimal perturbations, contain streamwise vortices. Tran- 
sient growth is closely associated with the sensitivity of the spectrum to small perturbations of 
the operator, quantified by the pseudospectra [18]. 

The concepts of pseudospectra, non-normality, transient growth, and optimal perturbations, 



have been rarely applied to Taylor-Couette flow [21, |22|] , probably because conventional linear 
stability analysis is so successful in explaining the transitions observed. However, the large 
parameter space of Taylor-Couette flow offers the possibility of tuning the system from a nor- 
mal operator, which does not support any transient growth, to a highly non-normal operator. 
Because of translational and Galilean invariance, plane Couette flow depends only on the dis- 
tance 2d between the plates, the relative velocity 2AU of the plates, and the kinematic viscosity 
v, which are combined into a single nondimensional parameter, the Reynolds number defined 
conventionally as Re = AUd/u. In contrast, because of curvature effects in Taylor-Couette 
flow, both the inner and outer radii r m and r out (or gapwidth 2d = r out — r- m and average 
radius f = (r; n + r out )/2) play a role. Because transformation to a rotating reference frame 
would introduce a Coriolis term into the equations, the flow depends on the angular velocities 
of both the inner and the outer cylinders, f2j n and f2 ut- The five dimensional parameters can 
be combined into three non-dimensional parameters in various ways. Our choice is rj = r; n /r out , 
fx = r2 ou t/r2i n , and Re = r ir Sl- m d/v. 

The idea of trying to approach the stability of plane Couette flow via Taylor-Couette flow is 
an appealing one and has inspired a number of investigations. The axisymmetric Taylor-vortex 
solution undergoes a secondary bifurcation to a non-axisymmetric wavy- vortex solution [Q, |6| . In 
a search for states intermediate in complexity between the Couette solution and turbulent plane 



Couette flow, Nagata [23, 24] took the limits of a narrow gap and almost corotating cylinders 
and discovered that, while the Taylor-vortex solution ceases to exist as the Coriolis term (i.e. 
the average angular velocity) is decreased, the wavy-vortex solution could be continued to the 
plane Couette limit. Faisst & Eckhardt |25| showed that these wavy solutions also exist for 
counter-rotating cylinders and that the lowest Reynolds number at which they first appear (via 
a saddle- node bifurcation) becomes independent of the rotation ratio as rj — ► 1. Finally, Prigent 
and Dauchot [Eq] have discovered that an analog of the spiral turbulence state of Taylor-Couette 
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flow also exists for plane Couette flow. 

Our investigation goes somewhat in the opposite direction. While the investigations cited 
above considered steady or traveling states known to exist in Taylor-Couette flow and continued 
them to the plane Couette limit, we take the ideas of pseudospectra, non- normality, transient 
growth, and optimal perturbations, apply them to Taylor-Couette flow and investigate the limit 
as Taylor-Couette flow approaches plane Couette flow. We focus on a given representative 
azimuthal and axial wavenumber and explore the stability characteristics of counter-rotating 
Taylor-Couette flow for a range of Reynolds numbers and radius ratios. In a related study, 
Meseguer |2^] studied transient effects - optimized over azimuthal and axial wavenumbers - 
for a specified radius ratio and varying Reynolds number and angular velocity ratio. Since 
Taylor-Couette flow is linearly stable for f2j n = 0, Meseguer suggests that transient growth 
plays an important role in the transition to turbulence observed experimentally pi for f2; n = 0, 



— O out 3> 1, a theory contested by Gebhardt & Grossman [21]. 

Why calculate transient growth in Taylor-Couette flow, one of the major success stories of 
linear stability theory? Transient growth is a new tool in linear stability theory, whose worth 
and longevity have yet to be proved. It has up to now been applied to flows in which transition 
from the basic laminar state is not understood. It is therefore worthwhile to calculate transient 
growth in Taylor-Couette flow, in which transition is believed to be completely understood by 
hydrodynamicists. This could provide information about the value and significance of transient 
growth calculations. The other side of the coin is that the paradigmatic status of Taylor-Couette 
flow should be maintained. As new tools or measurements become available, they should be 
brought to bear on this basic flow in order for hydrodynamicists to extract additional knowledge 
about Taylor-Couette flow. 
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II. GOVERNING EQUATIONS AND NUMERICAL METHODS 

A. Taylor-Couette flow 

The Couette solution is the unique solution of the form U_ = Uc{r)eg to the incompressible 
Navier-Stokes equations: 

?= + (U-V)U= -VP+^-AU, (la) 
at Re 

V-U = 0, (lb) 

U(r = f-l) = l, U( r = f+l) = ^. (lc) 

V 

Distances have been nondimensionalized by (r out — r- m )/2 and velocities by r- m £l m . We recall 
that 

i] = r in /r out , fi = fi out /Oi n , Re = r in fL in (r out - r in )/ (2u) (2) 

where r; n , r out , Oi n , S7 ut are the inner and outer cylindrical radii and angular velocities, respec- 
tively, and v is the kinematic viscosity. Definition (|2|) of Re, chosen for compatibility with the 
plane Couette convention for the case fj, = —1, differs by a factor of two from the conventional 
definition employed in Taylor-Couette flow. The average radius is f = (1 +rj)/(l — rj). Thus, 
the limit rj — > 1 corresponds to f —* oo; we refer to either of these limits as convenient. 

The Couette solution is: 

U c (r) = Ar + -, (3a) 
r 

A= , B=, 2 ^r M > (3b) 

277(1+7/)' (l-77)(l- 77 2 ) 1 ' 

The expressions for A and B again differ from the standard ones by factors of two due to our 
use of (r ou t — Hn)/2 as unit of length. Expression (|3|) can be viewed as a superposition of solid 
body rotation, Ar, and the flow due to a point vortex, B/r. As rj — ► 1, the leading terms of 
these two contributions become equal and opposite: Ar ~ —B/r ~ (/j — l)f/4. For rj ~ 1, we 
therefore rewrite (||) in the following equivalent form not subject to this cancellation error: 

U c (r) = [ 2 ^ 3 ((/* + 1) + (M - 1)?/) + r 2 (3(/x - 1) + 4(/x + l)y + (/x - l)y 2 ) 

+ 2fj/ (Ou - 1) + (/x + l)y) - (/j - 1)(1 - y 2 )] (4) 



where 



(5) 
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The equations we consider are the Navier-Stokes equations (||) linearized about the Couette 
solution d) or (|): 

8u T | U c du r 2U C u _ )= dp | 1 _^_]l^£\ (6 a ) 

8t r 89 r 8r Re\ r r 2 r 2 89 / ' 

due_ + C/c5ug + Uc_ u + jj c 'u = --— + —( Aug - — + ——) (6b) 
<9i r 89 r 1 r rd9Re\ r 2 r 2 89 ) ' 

^ + VW = -ai + ^( A N' (6c) 

^ + -, + -^ + ^ = 0, (6d) 
subject to the boundary conditions 

u T = ug = u z = at r = f ± 1 (6e) 

with 

_ <9 2 19 1 9 2 a 2 



The Taylor-Couette geometry and the Couette solution are homogeneous in the azimuthal (9) 
and axial (z) direction, which are analogous to the streamwise (x) and spanwise (z) directions in 



plane Couette flow. Table [II A| summarizes the equivalence of the coordinate systems for plane 
Couette and Taylor-Couette flow. Solutions to (|6|) which are spatially bounded are therefore 
trigonometric in each of these directions, with wavenumbers m and f3, 

(u r ,u e ,u z ,p) = (u r ,u$,u z ,p) (r, t) exp (im9 + i[3z) . (6) 



The parameter space is too vast to permit full exploration. In this study, we limit ourselves 
to [i = —1, m = 0, and (3 = ir/2. We choose /i = —1 so that the average angular velocity 
vanishes. 

The choice m = of axisymmetric perturbations is made for simplicity. In plane Couette flow, 
the transient growth achieved by streamwise-independent perturbations is, while not maximal, 
very close to the optimal value. Although the study of transient growth in Taylor-Couette flow is 
relatively unexplored, its linear instability has been extensively studied, e.g. [|], ||, [?], |j| The 
linear instability undergone by Taylor-Couette flow with counter-rotating cylinders is usually 
non-axisymmetric and leads to spirals rather than vortices. However, the threshold associated 
with axisymmetric perturbations is very close to the actual non-axisymmetric threshold. Indeed, 
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Plane Couette 


Taylor- Couette 


X 


streamwise 


a 





azimuthal 


m/r 


z 


spanwise 


(3 


z 


axial 




y 


normal 




r 


radial 





TABLE I: Equivalence of coordinate systems for plane Couette and Taylor-Couette flow. 



the thresholds are so close that the instability was thought to be axisymmetric until 1966, when 
calculations by Krueger, Gross and DiPrima [||] confirmed experimentally by Coles [0] showed 
the first instability to be non-axisymmetric for \x sufficiently negative, more precisely [i < 0.78 
in the narrow-gap limit. Note that the correspondence between the streamwise wavenumber a 
of plane Couette flow and the azimuthal wavenumber m of Taylor-Couette flow is m ~ fa, since 
x ~ r9. Setting m to correspond to a fixed non-zero value of a would thus require increasing m 
through integer values as rj — > 1 or f — > oo. 

The choice P = ir/2 corresponds to an axial wavelength of 4. Since our choice of length scales 
dictates that the radial gap is of width 2, this means that a single vortex has the same axial 
as radial extent, i.e. is approximately circular. The axial wavelength corresponding to linear 
instability is somewhat smaller than this for counter-rotating cylinders by as much as 20%, but 
is nonetheless close. 



We eliminate u z and p by using (|6q ) and the condition of incompressibility (pelf), obtaining 
evolution equations in u = (u r ,ug). The perturbations (u r ,ug) are represented as series of 
Chebyshev polynomials in y = r — f, which are evaluated at the Gauss-Lobatto points 

The evolution equation is written symbolically for the state vector u = (u r ,ug) T as 

^-u = —iCu. (7) 
dt~ w 

The axial velocity component u z is calculated from (u r ,itQ) via incompressibility. 

Our codes for discretizing the Taylor-Couette operator were based on a pseudospectral rep- 
resentation of the radial and azimuthal velocity in the inhomogeneous normal direction, similar 



to the Matlab code for plane Couette flow written by Reddy [19, 28 1 and published in [29]. We 
tested our Taylor-Couette code in a number of ways. First, we took the limit 77 — > 1 while main- 
taining a fixed Reynolds number and verified that we obtained the same asymptotic growth rates 
as in plane Couette flow. We then took 77 — > 1 and verified that the critical Reynolds numbers 
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were those obtained by Krueger et al. ||]. Finally, for arbitrary values of 7], we verified that the 
critical Reynolds numbers were those given in Q. 

B. Optimal Growth and Pseudospectra 

As a measure of growth, we use the energy norm defined by: 



E(u) = \\u\\ E = I I \u r \ + — ^\Vu r \ + -y2~^\ue\ rdr+ \u \ r dr . (8) 

r=f-l V P P r J Jr=f-1 



E r ,z Eg 

with T>f = ^-§p{ r f)- The energy norm is particularly significant in hydrodynamics, because the 
nonlinear terms of the Navier-Stokes equations conserve energy. This lends significance to the 
study of the linearized equations (^), since any energy growth that takes place must occur via 
linear mechanisms [29]. The maximal energy growth at time t for u evolving according to ([?]) is 



defined by: 

«(o)^o IK0)||| 

Thus the maximal energy growth is given by the energy norm of the operator exp(— iCt). 

The norm of a normal operator is its dominant eigenvalue A max ; this is the largest factor 
by which matrix multiplication can increase the norm of a vector. For a non-normal operator, 
cross-terms between non-orthogonal eigenvectors typically contribute to the norm of a vector. 
Instead, it is the singular vectors which are orthogonal; the norm of the operator is given by 
the largest singular value <7 max . The singular value <J max and its corresponding normalized right 
and left singular vectors « max , w m ax satisfy: 

exp(-i£t)u max = o- max u max (10) 

i.e. linear evolution from initial condition u max leads to the state cr max ?; max . 

The statements above are all inner-product dependent: an operator is normal or not and 
vectors are orthogonal or not with respect to a particular inner product. The singular value 
decomposition is inner-product dependent as well. Since we investigate growth in the energy 
norm, we seek the largest singular value and its corresponding left and right singular vectors in 
the energy norm as well. Standard software, however, provides singular value decompositions 
with respect to the 2-norm. Additionally, each of the values on the Gauss-Lobatto grid points 
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must be multiplied by a weight appropriate for calculating the energy (|8|). We compute the 
elements My of an N x N Hermitian matrix M by taking the inner products, derived from (||), 
between two eigenfunctions <3?j and of C The Cholesky factorization of M = F^F is then 
used to convert the energy norm of the operator exponential to an L2-norm of the weighted 
matrix exponential according to 



exp(-i£t)||| « ||Fexp(-iDt)F _1 || 2 (11) 



with D as an N x N diagonal matrix consisting of the eigenvalues of C (see [ 29 1 for more details) . 



We calculate the largest singular value (under the 2-norm) of the operator on the right-hand-side 



of (11). Its square is the maximal energy growth. The number N of eigenfunctions has been 



chosen large enough to ensure converged results. 
The optimal growth is defined as 

G max =supG(t). (12) 

t>o 

If C has an eigenvalue with positive imaginary part, then ||u|| grows exponentially in time (for 
any norm) and so G max = oo. Thus calculations of optimal growth are meaningful only for 
operators which are linearly stable. 

We also wish to keep track of the quantities responsible for achieving the maxima in (^) and 
(12). The time for optimal growth £ max is that which achieves the maximum (sup) in (|i~2|), i.e. 



G max = G(t max ). The optimal input, denoted by u(0), is the normalized initial condition which 
achieves the maximum (sup) in (||) for t = i max . The optimal output, denoted by u(t max ), is 
the velocity field resulting from the linearized Taylor-Couette evolution ([/]) starting from the 
unit-energy optimal input u(0); its energy gain is G max . 



A non-normal operator C is also characterized by its pseudospectra Pq| . The e- 
pseudospectrum A £ (£) is the set of complex values z (parametrized by e) which satisfies the 
property: 



KzI-C^W > e- 1 



(13) 



Equivalent definitions of the e-pseudospectrum are given in IS, 19, 28, 29, 3C]. The definition 
of the pseudospectra, like that of transient growth, depends on the norm or inner product. For 
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a normal operator, the e-pseudospectrum is the union of the balls of radius e surrounding each 
eigenvalue. For a non-normal operator, on the other hand, the e-pseudospectrum may be much 
larger. 

Kreiss ' theorem relates the optimal growth and the pseudospectra by the following inequality: 



Gmax > sup e 1 sup lm(z) . (14) 
oo V zeA e (£) I 



The right-hand-side of (14) maximizes (over all strictly positive values of e) the ratio of the 
distance to the real axis of any point in the e-pseudospectrum A t to the value of e. In practice, it 
is found, both for Taylor-Couette flow and for plane channel flows, that the optimal growth given 
by the left-hand-side of (|l4|) is approximately twice the lower bound given by the right-hand-side 
of ([1|). 

For computing the optimal growth, we used the Matlab code written by Reddy |l^, |2^| , given 
in 1 29]. For computing the pseudospectra, we used the code EigTool written by Wright |31], 



which in turn makes use of the algorithm developed by Trefethen [30] and is available at website 



tittp : //www . comlab . ox . ac . uk/pro j ect s/pseudospectra 



III. RESULTS 



We begin by presenting the optimal growth as a function of r/ and Re. We recall that 
throughout our study, we fix the azimuthal wavenumber m = a = 0, the axial wavenumber 
P = 7r/2, and the angular velocity ratio fi = —1. Figure |l] shows the contours of constant 
optimal growth G max . Inside the shaded region, Couette flow is linearly unstable to axisymmetric 
perturbations. The boundary of this region is the critical Reynolds number ReL(rj). Rei — ► oo 
as 77 — ► 1, as expected since plane Couette flow is linearly stable for all Reynolds numbers. We 
may also consider the rightmost portion of the linear stability boundary as a function r/ii^Re). 

A striking feature is that the maximum growth for a fixed Reynolds number, indicated by 
the triangles in figure [l], is always achieved for a radius ratio r] = rj op t(Re) which is less than 
one. We propose a possible explanation for this trend. As we decrease the radius ratio 77, the 
asymptotic growth rate, i.e. the imaginary part of the least stable eigenvalue, increases. At 
the same time, the non-normality of the operator decreases, resulting in diminished transient 
growth. The combination of these effects results in a maximum growth rate that is achieved 
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FIG. 1: Contours of optimal growth for Taylor-Couette flow in the (rj,Re) plane. Shaded area indicates 
the region of linear instability. Triangles indicate r] pt(Re), the value of r\ at which maximum growth is 
attained for a given value of Re. 

for values of rj less than one. We see that r/ op t ~ 0.6 for Re = 20, increases to a maximum of 
T/opt ~ 0.96 for Re ~ 110, and then abruptly decreases and terminates by meeting the linear 
instability boundary r/i ~ 0.9 at Re ~ 130. For Re > 130, the maximum growth is achieved for 
rj = rjL(Re). The enlargements in figure ^ show the typical behavior of the contours of G max 
for Re near 130 and for Re near 300. For Re < 130, the optimal growth is fairly weak, varying 
from a factor of 1 to 21. Arbitrarily high values of G max can be attained by increasing Re, since 
for plane Couette flow, i.e. rj = 1, it is known [jl!| that G max ~ Re 2 . In fact, over the range 
300 < Re < 310, G max is approximately 20% higher for r\ = rj^ than for rj = 1. In this range, 
the G max contours are nearly vertical as they approach the linear instability boundary, meaning 
that G max is far more sensitive to a decrease in rj than to an increase in Re. It is near Re = 310 
that plane Couette flow undergoes a sudden unexplained transition to turbulence. Table |Tj gives 
selected numerical values of G max . 

We now study in detail two contrasting cases: rj = 0.5, Re = 125 and 7] = 0.99, Re = 350. 
In the first case, rj = 0.5 or equivalently r = 3, curvature obviously plays an important role. 
The second case, rj = 0.99 or equivalently f = 199, is very near the plane Couette limit. The 
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FIG. 2: Close-up views of contours of optimal growth. Above: for 120 < Re < 130, each contour changes 
slope at a value rj op t(Re) an 0.96. Below: for Re > 130, G max increases as rj decreases to the linear 
stability boundary rj^. 

Reynolds numbers have been chosen to be close to the linear instability threshold Rei in each 
case in order to maximize transient growth while remaining within the linearly stable region. 

Figures |3| and |I| show the optimal input u(0) and output ;u(i max ) for each case. The least 
stable eigenvector, i.e. that with the smallest decay rate, is not shown, but resembles the 
optimal output. The upper portion of each figure shows the meridional velocity fields (u r ,u z ) 
of the optimal input and output, while the lower portion shows contours of azimuthal velocity 
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V 

Re 


0.5 


0.6 


0.7 


0.8 


0.9 


0.95 


0.99 


0.999 


0.9999 


Plane 
Couette 


50 


1.95 


2.39 


2.85 


3.23 


3.38 


3.34 


3.26 


3.23 


3.23 


3.23 


75 


2.62 


3.51 


4.66 


5.95 


6.90 


7.02 


6.90 


6.85 


6.84 


6.84 


100 


3.29 


4.69 


6.74 


9.40 


11.70 


12.14 


12.01 


11.93 


11.92 


11.92 


125 


4.17 








18.45 


18.87 


18.59 


18.46 


18.44 


18.44 


150 


5.53 










27.72 


26.65 


26.44 


26.42 


26.42 


300 














111.60 


104.87 


104.74 


104.73 


500 
















291.92 


290.39 


290.35 



TABLE II: Optimal growth G max in the narrow-gap limit. Dashed entries indicate linear instability, i.e. 
G max = oo. Note that, for fixed Re, the maximum optimal growth is achieved for rj < 1. 



ug. The meridional velocity fields consist of vortices whose axes are oriented in the azimuthal 
direction, similar to the eigenvectors which lead to Taylor vortices at slightly higher Reynolds 
numbers and to the streamwise vortices which are the optimal inputs in plane Couette flow. The 
azimuthal components of the optimal inputs and outputs, shown in the lower portions of figures 
^ and are in phase opposition with the vortices, with nodal lines at z = 1,3 going through 
the vortex centers. 

Figure || shows the evolution in time of the energies in the meridional components E r>g and 
in the azimuthal component Eg starting from the optimal input u(0). While the optimal inputs 
u(0) are concentrated primarily in the meridional components, it is the azimuthal component 
which dominates the optimal outputs u(t max ). For this reason, in order to show the qualitative 
geometric features of the two fields, the inputs and outputs of Figures || and || use different scales 
for the arrow lengths and for the contour levels. This evolution corresponds to the generation of 
streaks - deformations of the azimuthal velocity profile - by the vortices. This physical process, 



referred to as the lift-up mechanism, has been described, e.g., in [32] and is believed to be a key 
element in the transition to turbulence in plane Couette flow. 

The graph on the left of figure |5| shows the evolution of Eg and E r<z starting from the optimal 
input u(0) and from the least stable eigenvector for r] = 0.50, Re = 125. The initial energies are 
Eg(0) = 0.2 and .EV iZ (0) = 0.8. Initially, over < t < i max 5, Eg rises while E r>z decreases, 
attaining values of Eg(t max ) m 4 and E r ^ z (t maiX ) 0.1, with a ratio Eg / E r>z (t max ) ~ 40. Over 
the interval 5 < t < 8, E r z increases, while Eg continues to decrease. With further evolution, 
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FIG. 3: Optimal perturbations for r\ = 0.50, Re = 125. Left: optimal input. Right: optimal output. 
Above: Meridional velocity field (u r ,u z ). Below: contours of ug. The energy of the input (output) is 
primarily in the meridional (azimuthal) components. For this reason, arrow lengths and contour levels 
are scaled differently for the input and the output. 

both energies decrease as u(t) converges towards the least stable eigenvector. From their values 
at t = 20, we estimate E 6 /E rjZ (oo) ki 0.73/0.023 ki 32. 

The energy evolution for the case r\ = 0.99, Re = 350, shown on the right of figure [|, resembles 
that for plane Couette flow. The optimal input is almost exclusively meridional, with negligible 
azimuthal component: E rz ~ 0.996 while Eg(0) ~ 0.004. By i max ~ 66, Eg has increased to 164 
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FIG. 4: Optimal perturbations for r) = 0.99, Re = 350. Left: optimal input. Right: optimal output. 
Above: Meridional velocity field (u r , u z ). Below: contours of ug. The energy of the input (output) is pri- 
marily in the meridional (azimuthal) components; arrow lengths and contour levels are scaled differently 
for the input and the output. 

and E r:Z decreased to 0.14, a ratio Eg / E rtZ (t mSLX ) of 1170 for the optimal output. This ratio is 
approximately maintained as both energies slowly decrease during the evolution of u(t) towards 
the least stable eigenvector. 

In the case r\ = 0.5, Re = 125, two arrays of vortices are present in the optimal input, a larger 
and stronger array near the inner cylinder and a smaller and weaker array near the outer cylinder. 
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FIG. 5: Time evolution of optimal growth. Solid curves denote azimuthal energy Eg(t), long-dashed 
curves meridional energy E riZ (t) during evolution from optimal input u(0). Higher and lower short- 
dashed curves represent Ee{t) and E r ^ z {t), respectively, for least stable eigenvector. Triangle corresponds 
to tmax- Above: r\ = 0.50, Re = 125. Below: rj = 0.99, Re = 350. 
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FIG. 6: Square of angular momentum of basic Couette solution (rUc) 2 - Solid curve: rj = 0.5, axis on 
left. Dashed curve: rj — 0.99, axis on right. Rayleigh's criterion for inviscid instability states that Uc is 
unstable where (rUc) 2 is a decreasing function of r. 

Three arrays are present in the optimal output, whose radial extent and strength decreases in 
going from the inner to the outer cylinder. In the case r\ = 0.99, Re = 350, the optimal input 
contains one array of vortices and the optimal output contains a second weaker, narrower array 
near the outer cylinder. Some light can be shed on the form of these perturbation fields and 
on the difference between the two cases by Rayleigh's criterion for instability in Taylor-Couette 
flow. 

Rayleigh's argument, valid for inviscid and axisymmetric flow, is that perturbations inter- 
changing rings of fluid at different radii (e.g. Taylor vortices) will be favored or opposed by the 
ambient pressure gradient, according to whether the square of the angular momentum (rUc) 2 
decreases or increases radially outwards. For counter-rotating cylinders {ji < 0), the sign of 
d(rUc) 2 /dr changes within the gap. Rayleigh's criterion is then applied to argue that only the 
inner portion of the gap is unstable. This modification is justified by three related tendencies 
H . First, the unstable eigenvector is concentrated near the inner cylinder, where d(rUc) 2 / dr is 
negative. Second, the axial wavelength corresponding to the most unstable or least stable eigen- 
vector decreases, favoring vortices which remain closer to circular. Third, the critical Reynolds 
number for linear instability increases, meaning that the critical Reynolds number based on the 
unstable portion of the gap remains nearly constant. 

Figure ^ shows that the square of the angular momentum decreases radially outwards over 
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FIG. 7: Pscudospcctra for Taylor-Couette flow. Above: i] = 0.5 and Re — 125. Contours corre- 
spond to e = lO" 18 , 10" 15 , . . . , 10°. Below: 77 = 0.99 and Re = 350. Contours correspond to 
e = 10- 2 - 4 ,10- 2 -\...,10°. 

the interval r < f — 0.5 for 77 = 0.50 and over the interval r < f for 77 = 0.99. Although exact 
application of Rayleigh's criterion would lead to optimal perturbations far more concentrated 
near the inner cylinder than they actually are, the criterion provides a heuristic explanation 
for the asymmetry. Rayleigh's criterion is usually invoked to explain linear instability, i.e. 
exponential growth. However, a modified version of the criterion should apply to transient 
growth as well. 
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Finally, we show the spectrum and pseudospectra of the operators C for the two cases in 
figure |]. Both pseudospectra plots contain contours (for the smaller values of e) surrounding 
individual eigenvalues and bulb-shaped contours (for the larger values of e) surrounding the 
entire spectrum. For both values of r/, the contours for a fixed small value of e surrounding the 
eigenvalues near the real axis are wider than the contours surrounding the eigenvalues farther 
from the real axis. The bulb-shaped contours, however, differentiate between the two values of 
rj. The spectrum for rj = 0.5 contains eigenvalues near the real axis with real parts extending 
to approximately ±0.7; a bulb-shaped pseudospectral contour thus remains a fairly constant 
distance from the spectrum. The spectrum for r/ = 0.99 is, in contrast, quite localized on the 
imaginary axis; in this case, a bulb-shaped pseudospectral contour protruding into the unstable 
half-plane is a first indication of non- normal effects. 

We use more detailed calculations of the pseudospectra to compute approximations to the 
lower bound (|i~4| ) on optimal growth of Kreiss' theorem. For r\ = 0.5, we obtain an upper 
bound Im(z) < 0.1 for the e-pseudospectrum with e = 10~ L15 , which yields the lower bound 
G max > (10 115 xO.l) 2 = 1.99, about half of the exact value G max — 4.17 that we have calculated. 
For rj = 0.99, we obtain an upper bound Im(z) < 0.03 for the e-pseudospectrum with e = 10~ 2 ' 45 , 
which yields the lower bound G max > (10 2 - 45 x 0.03) 2 = 71.49, again about half of the exact 
value of 155. 

IV. CONCLUSIONS 

We have calculated the pseudospectra and optimal transient growth for Taylor-Couette flow. 

Our major result is that, for a fixed Reynolds number, the optimal transient growth is 
achieved for a radius ratio r/ = rj opt < 1, rather than for the plane Couette limit rj = 1. This is 
due to the combined effect of increasing modal growth and decreasing non-normality. As shown 
in figure EL for Re < 130, the optimal transient growth for fixed Re occurs at an intermediate 
value of T) ranging from 7] opt = 0.6 at Re = 20 to 77 op t = 0.96 at Re = 110. For Re > 130, the 
optimal transient growth for fixed Re increases as r\ is decreased from the plane Couette limit 
r] = 1 to the linear stability boundary rj = 7]l, as shown in figure ^. 

If transient growth on the order of G max ~ 100 is suspected to initiate transition in plane 
Couette flow near Re = 300, then transition should also occur near Re = 300 outside the linear 
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instability domain in Taylor-Couette flow, albeit in very narrow gaps {rj > 0.986). If these 
predictions fail to hold, then theories concerning transition initiated by transient growth must 
be modified or even abandoned. Another possible direction for future studies concerns regions 
of (rj, Re) in which transient growth competes with linear instability; such situations should be 
included in a comprehensive theory of transition initiated by non-normal effects. 

Two cases have been considered in detail: rj = 0.99, Re = 350 and 77 = 0.5, Re = 125. Both 
show transient growth. As could be expected, the case r\ = 0.99, Re = 350, shows much higher 
transient growth, both because of its greater resemblance to plane Couette flow and also because 
of its higher Reynolds number. As shown in figure |5|, for r\ = 0.99, Re = 350, energy grows from 
the optimal input by a factor of G max = 164, while for rj = 0.5, Re = 125, we have G max = 4. 

As is the case for plane Couette flow, the physical mechanism accompanying transient growth 
consists of the conversion of vortices (here, termed azimuthal rather than streamwise) into 
streaks, i.e. perturbations of the basic Couette profile. This is seen in the ratio of azimuthal 
to meridional energy in figure ||| of the optimal inputs u(0) and outputs u(t max ). The optimal 
outputs, depicted in figures ^ and |4|, are of higher amplitude near the inner cylinder, especially 
in the case rj = 0.5, Re = 125, as indicated by Rayleigh's criterion for centrifugal instability 
illustrated in figure |6| 

Although both sets of pseudospectra in figure ^ show signs of non-normality, those for rj = 
0.99, Re = 350 show more deviation from the spectrum than those for rj = 0.5, Re = 125. 

Our study is restricted to the azimuthal wavenumber m = a = 0, axial wavenumber (3 = n/2 
and angular velocity ratio ji = — 1. Previous results on Couette flows lead us to believe that 
these values of (a, (3) may be representative of a fairly large portion of parameter space, since 
both the maximal optimal growth rate in plane Couette flow and the maximal linear growth 
rate in Taylor Couette flow are close to those achieved for a = 0, [3 = tt/2. Varying ji, however, 
leads to major qualitative changes, as it does in other aspects of Taylor-Couette flow. Studies 
encompassing a wide range of parameter values [p2| are clearly desirable. 



In the words of Faisst & Eckhardt [25], Taylor-Couette flow provides an "embedding" of 
plane Couette flow. Numerical simulations and laboratory experiments have extensively docu- 
mented the way in which plane Couette flow changes as its single non-dimensional parameter, 
the Reynolds number, is increased. Taylor-Couette flow provides an ensemble of other param- 
eter paths along which to approach or to step back from plane Couette flow. Our hope is 
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that this preliminary study and that of [22| of transient growth and pseudospectra in Taylor- 
Couette flow, will help to increase understanding of both Taylor-Couette flow and of the effects 
of non-normality. 
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